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ABSTRACT 

We report a dual-band observation at 223 and 654 GHz (460 /xm) toward 
an ultracompact (UC) Hll region, G240.31+0.07, with the Submillimeter Array. 
With a beam size of V!5 x 0^/8, the dust continuum emission is resolved into two 
clumps, with clump A coincident well with an H2O maser and the UC Hll region. 
The newly discovered clump, B, about V!3 {c^ 8.3 x 10^ AU) to the southwest 
of clump A, is also associated with H2O masers and may be a more recent star- 
forming site. The continuum flux densities imply an opacity spectral index of 
/? = 1.5 ± 0.3, suggestive of a value lower than the canonical 2.0 found in the 
interstellar medium and in cold, massive cores. The presence of hot {c^ 100 K) 
molecular gas is derived by the brightness ratio of two H2CO lines in the 223 GHz 
band. A radial velocity difference of 2.5 ± 0.4 kms~^ is found between the two 
clumps in C^^O (6 — 5) emission. The total (nebular and stellar) mass of roughly 
58 Mq in the central region is close to, but not by far larger than, the minimum 
mass required for the two clumps to be gravitationally bound for binary rotation. 
Our continuum data do not suggest a large amount of matter associated with the 
H2 knots that were previously proposed to arise from a massive disk or envelope. 

Subject headings: circumstellar matter — stars: individual (G240.31+0.07) — 
stars: pre-main sequence — stars: early-type — submillimeter 

1. INTRODUCTION 

The bright far-infrared source, G240.31+0.07 (hereafter G240), at a distance of 6.4 kpc 
(MacLeod et al. 1998), is known to harbor an ultracompact (UC) Hll region associated 
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with H2O and OH masers (Migenes et al. 1999; Caswell 2003). The flux density at 100 /xm 
suggests a luminosity of L = 10"^'^ Lq^ consistent with a spectral type 08.5 zero-age main 
sequence (ZAMS) star (MacLeod et al. 1998). Observations with the Very Large Array 
(VLA) detected a point source with radio free-free emission corresponding to ionization 
by a ZAMS BO. 4 star (MacLeod et al. 1998). On the large scale, the CO (1-0) emission 
exhibits broad line wings up to 38.6 kms~^, possibly tracing a bipolar outflow (Shepherd 
& Churchwell 1996). With improved angular resolution of 20'^, Hunter (1997) resolved the 
CO (3-2) emission and found a prominent bipolar outflow at a position angle (P. A.) of 138° 
and a possible second outflow axis at P. A. 101°. Each of the two bipolar outflow extends 
roughly 80^' in length. Subsequent infrared study found two bright, elongated H2 emission 
features within 16^' (0.5 pc) of the UC Hll region (Fig. [U Kumar et al. 2002). Kumar 
et al. (2003) analyzed emission lines of H2 and [Fell] and concluded the H2 emission being a 
result of non-dissociative J-shocks. The C^^O (2 — 1) emission, from 65 to 70 kms~^, shows 
a linear velocity gradient of 0.1 kms~^ arcsec"^ along the H2 axis with a 33^' beam. The 
authors further argued that the H2 emission features arise in a rotating disk or envelope 
around G240. 

Using single-dish telescopes, millimeter and submillimeter observations toward G240 
can not resolve structures smaller than 12'^ Interferometric observations with high angular 
resolutions are crucial to investigating the linear velocity gradient seen in the C^^O (2 — 1) 
emission. Our current study is enabled by the Submillimeter Array^ (SMA) (Ho, Moran, & 
Lo 2004) with its highest frequency band at 654 GHz. Dust continuum emission and the 
higher transition of C^^O (6 — 5) {E^p = 110.6 K) should better trace warm components in 
this region. Moreover, the temperature and the dust properties can be studied by combining 
the simultaneous observation in the 223 GHz band. 



2. OBSERVATIONS AND DATA REDUCTION 

The observation was carried out on Feburary 16, 2005, using 6 elements of the SMA in 
the dual-band mode, which allows simultaneous operation of the 223 and 654 GHz receivers. 
The phase center was toward (7^44"^52^0, — 24°7'42'.'6)(J2000). The observing cycle comprised 
scans of 0736+017, Titan, VY CMa, and G240, and was repeated every 25 minutes. The 
total on-source integration time was about 80 minutes for G240. Data inspection, bandpass 
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and flux calibrations were done within the IDL superset MIR. The Jovian moon, CaUisto, 
was used to cahbrate the passband as weU as the flux density scale. The flux measurements 
should be accurate to within 30% and 15% at 654 and 223 GHz, respectively. Details of the 
observations are listed in Table [H 

A major difliculty for interferometric imaging at 654 GHz is the lack of sufliciently 
strong quasars that can be used to derive complex gains. Our strategy is to utilize the 
nearby (5° away from G240) H2O masers at 658.00655 GHz (Menten & Young 1995) in 
VY CMa to trace atmospheric phase fluctuations. The H2O masers were observed in the 
same sideband with the C^^O (6 — 5) transition. Calibration of the 654 GHz band was done 
in two steps. First, phase-only gain solutions were derived from the H2O masers, whose flux 
density, however, may change on the timescale of a few hours due to its possibly polarized 
and highly variable nature. Without a stable flux density, the masers can not be used for 
gain amplitude calibration. Instead, we examined the continuum emission of VY CMa and 
found it to be a structureless source of 7 Jy with a 1^.'5 x O^.'S beam. The peak position of 
the continuum emission agrees with the maser spot to within 0^.^06, signiflcantly smaller than 
the beam. We therefore derived complex gains from the continuum emission of VY CMa to 
correct temporal amplitude variation. Gain solutions from the two steps were then multiplied 
together and applied to the other sources. Judging from the peak poisitons of Titan (40° 
away) and 0736+017 (20° away), we estimated the astrometry should be good to within 0'.'3. 
Both sidebands, each of a 2 GHz bandwidth, were used to generate a line-free continuum 
map in a mult i- frequency synthesis, resulting a frequency equal to the average of the upper 
and lower sideband frequencies. Imaging was performed with the MIRIAD package. 

Calibration of the 223 GHz band followed the conventional process by solving complex 
gains against the quasar 0736+017, whose flux density was assumed to be 1.9 Jy. The gain 
solutions were then transferred to the other sources. A line- free continuum map and channel 
maps of G240 were generated with uniform weighting. 

3. RESULTS AND DISCUSSION 

3.1. Morphology and Kinematics 

The 654 GHz continuum emission of G240 is resolved into two clumps (Fig. [1]), with 
clump A coincident well with a VLA 6 cm point source (Hughes & MacLeod 1993) and 
an H2O maser. The newly discovered clump, B, is also associated with H2O masers and 
is probably a more recent star- forming site (Codella et al. 2004). The two sources may 
be responsible for the two bipolar outflows in CO (3 — 2) emission (Hunter 1997). The 
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peak of clump A is at (7^44^52^04 -24°0r42'/20)(J2000) with 2.76 Jybeam-i while that 
of clump B is at (7^44^51^96 -24°0r42^/92)(J2000) with 2.44 Jybeam-^ The projected 
angular separation between the two peaks is V!3^ equivalent to 8.3 x 10^ AU, with a P. A. of 
57°. No other continuum source is detected above the 3a sensitivity of 0.39 Jybeam"^ in 
the field of view. 

Fig. [2] shows the C^^O (6 — 5) integrated intensity map and the spectra toward clump A 
and B with Gaussian fits. Both clumps are associated with C^^O emission and show large 
Gaussian linewidths attributed to nonthermal eflFects, e.g. turbulent motions. The spectrum 
toward each clump has a Gaussian linewidth of roughly 1.7 kms~^, much larger than the 
thermal linewidth of 0.17 kms~^ for a clump at 100 K. Together with the geometric mean 
of 0^.^13 for the emission size, we estimate a virial mass of 68 Mq from the C^^O clump. The 
Gaussian peak velocity of clump A is va = 66.1 ±0.2kms~^, similar to the ambient velocity 
of 67.1 kms"-*^ in the CO (1 — 0) emission (Shepherd & Churchwell 1996). Clump B is at 
Vb = 63.6 ± 0.4 kms~^ based on our C^^O (6 — 5) spectrum and was not detected in the 
previous C^^O (2 — 1) observation with a larger 33^' beam (Kumar et al. 2003). Furthermore, 
our observation does not detect a velocity gradient along the east-west direction as what has 
been seen in the large-scale C-^^0 (2 — 1) emission. The kinematics at small scales (0.03 pc) 
do not follow the kinematics at large scales (1.2 pc). 



3.2. Line Emission in the 223 GHz band 

Formaldehyde (H2CO) possesses many transitions of greatly different excitation ener- 
gies and is a good probe for kinetic temperature (Mangum & Wootten 1993). Two H2CO 
transitions, 3o,3 — 2o,2 and 32,2 — 22,1, with upper state energy levels of Eup = 21 and 68 K, 
respectively, were observed in the 223 GHz band (Fig. [3]e and [3^). We calculate an average 
temperature in the central 5^' x 5^' region by comparing the mean integrated brightnesses 
of 13.8 ± 0.8 and 4.7 ± 0.8 Kkms~^ in the 3o,3 — 2o,2 and 32,2 — 22,1 lines, respectively. 
Assuming local thermodynamic equilibrium (LTE), the gas rotational temperature is about 
= 96 ± 35 K if both lines are optically thin. 

In addition to the two H2CO lines, other detected molecular lines are listed in Table [2] 
and velocity-integrated brightness maps are shown in Fig. [3l Assuming LTE and optically 
thin emission, we calculate the abundance for each species (Table [2j). Several lines, e.g. the 
two CN lines, are clearly affected by the lack of short baselines, suggesting that a significant 
fraction of the emission is diffuse. Although the low angular resolution at 223 GHz can not 
spatially resolve the two clumps, two distinct kinematic clumps are well separated in velocity 
space in the CH3OH and H2CO emissions. Fig. ^ shows the position-velocity diagram of 
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the H2CO emission as an example. 

The SiO {5 — A) v = emission does not appear to be associated with extensive outflow 
activities (Fig. [3]c). The derived abundance of 2.4 x 10~^^ is similar to the average ambient 
abundances of 10~^^ to 10~^^ in clouds with active star formation (Codella, Bachiller, & 
Reipurth 1999). There is no clear evidence for SiO abundance enhancement toward G240 
due to outflow activities. 



3.3. Core Properties and Mass Estimate 

The 223 GHz data indicate the presence of extended structure (i.e., a halo), as the 
visibility amplitudes increase with decreasing baseline length. To estimate the fluxes of the 
compact and extended components, we flt the visibilities with a "core-halo" model consisting 
of a point source and a one-dimensional Gaussian. This decomposition gives 223 GHz flux 
densities of 0.24 ± 0.04 and 0.37 ± 0.06 Jy for the core and halo, respectively, and a halo 
size of about Similar core-halo dust emission structures have been observed toward other 
massive star- forming regions (Su et al. 2004). The 654 GHz data is not sensitive to this 
spatially extended halo. 

If the dust emission is optically thin, the continuum flux density Fj^ = Biy{Td)/ D'^^ 

where /^^, M^^ Bj^iTd)^ and D are the dust opacity at the observing frequency, z/, the dust 
mass, the Planck function at a dust temperature, T^, and the distance to the source, re- 
spectively. The opacity spectral index (3 = dln/^^/rflnz/ can be measured by comparing flux 
densities at two frequencies. At 654 GHz, the Rayleigh- Jeans approximation produces 15% 
error at a temperature of 100 K and is not applied to our calculation. Because of different 
uv coverage at 223 and 654 GHz, we compare flux densities in maps made with visibilities of 
projected baselines within (22 — 52) kA, the range that the two frequencies have in common 
and corresponds to the compact component within the halo. The flux densities at 223 and 
654 GHz are 0.28 ± 0.04 and 10.4 ± 3.1 Jy, respectively, and imply that /? = 1.5 ± 0.3, 
smaller than the canonical value of 2 in the interstellar medium and in cold, massive cores 
(Hill et al. 2006). Small values of (3 have been observed toward several massive star- forming 
regions (Chen et al. 2006; Beuther et al. 2004) and may be attributed to grain growth or 
merely a feature of high-density environment (Chandler et al. 1995; Draine 2006). On the 
other hand, a small /? can also be a result of the opacity at 654 GHz becoming important. 
Assuming a dust opacity of /^654GHz = 10 (Anim/0.25)~^'^ = 4 cm^g~^ (Hildebrand 1983), 
the peak intensity of 2.76 Jybeam"^ at clump A corresponds to a dust column density 
0.03 gcm~^, translating to an optical depth tq^agrz — 0.13. Hence, the observed small /? is 
not due to large opacity at 654 GHz. 



-6- 



The dust mass can be calculated with = Fj^D'^ / tvjyBjy{Tg) = 0.47 Mq by assuming 
thermal equilibrium between dust and gas and a 223 GHz dust opacity /^223GHz = 0.8 cm^ g~^. 
Hence, the total nebular mass is about 47 Mq if the gas-to-dust ratio is 100 (Savage & Mathis 
1979). The estimated virial mass of 68 Mq (Sec. 13.11) is found to be slightly larger than the 
nebular mass. The cloud may be marginally self-gravitationally bounded or under additional 
pressure exerted by the outer halo. The free-free emission of a ZAMS BO. 4 star accounts 
for the minimum stellar mass of 11 Mq (Hanson, Howarth, & Conti 1997). Therefore, 
the enclosed mass in the central region should be about 58 Mq. Given the radial velocity 
difference Av = 2.5±0.4kms~^ and the projected separation r = 8.3 x 10^ AU, the minimum 
mass required to hold the two clumps in orbit about each other will be (At;)^ r/G = 60 Mq. 
The total (nebular and stellar) mass in the central region is close to, but not by far larger 
than, the minimum mass required for the two clumps to be in binary rotation. We can not 
rule out the possibility of the two clumps being members of a wide binary system. Studies 
of the spectral energy distribution suggested about 390 Mq of matter, mostly in the halo, 
associated with G240 (Kumar et al. 2003). The two clumps in the central region should be 
bound inside the gravitational potential well generated by the massive halo. 

Neither of our continuum maps shows elongated structures corresponding to the H2 
knots that have been proposed to trace a disk or envelope (Kumar et al. 2003). Assuming 
typical conditions P ^ 1 and T 40 K for a disk (Mannings & Sargent 1997), we use the 
223 GHz continuum sensitivity to place an upper limit on the visual extinction, Ay^ for this 
hypothesized disk. The 2a level of 4 mJy beam~^ corresponds to 0.8 Mq within a beam. The 
lack of continuum emission suggests that the visual extinction toward the H2 knots should 
be no larger than 4.8 mag. In contrast to a massive disk of Ay = 10.7 mag, our data do not 
suggest such a large amount of matter associated with the H2 knots. 

We are grateful to Ken Young, Rob Christensen, Alison Peck and the SMA team for 
making the 690 GHz campaign possible. We also thank Ed Sutton for the discussion on 
CH3OH transitions. 
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Table 1. Parameters for the Dual-band Observation 



Frequency band (GHz) 654 223 



Calibrator VY CMa 0736 + 017 

Primary beam (arcsec) 19 56 

Frequency coverage (GHz) 646.8 - 648.8 216.5 - 218.5 

656.9 - 658.9 226.5 - 228.5 

Typical T^ys (K) 2000 200 

uv range (kA) 22 - 150 8 - 52 

Largest visible scale (arcsec) 4 11 



Continuum 

Weight scheme 
Synthesized beam 
rms noise (mJybeam"^) 

Line 

Weight scheme Natural Uniform 

Synthesized beam 1'.'7 x I'.'O (0°) 5'.'0 x 2'.'7 (2°) 

rms noise (K) 5.3 0.2 

Channel width (km s"^) 0.6 1.2 



Uniform 
1'.'5 X O'.'S (6°) 
130 



Uniform 
4'.'4 X 2'.'4 (1°) 
2 
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Table 2. Detected Molecular Line Emission 



Frequency 


Species 


Transition 


E 


rppeak 
-'-B 








(GHz) 


(X) 




(K) 


(K) 


(Kkms-i) 


(cm~^) 




216.643303 


SO2 


222,20 ~ 22i^2i 


248 


0.8 


4.7 


1.5E+15 


4.0E-9 


217.104981 


SiO 


5-4t; = 


31 


1.2 


4.3 


8.8E+12 


2.4E-11 


217.238539 


DCN 


3-2 


21 


2.0 


6.0 


1.2E+13 


3.2E-11 


218.222186 


H2CO 


3o,3 — 2o,2 


21 


4.3 


23.5^= 


7.1E+14 


2.0E-9 


218.440050 


CH3OH 


42,2 ~ 3i,2 E 


45 


3.0 


10.1 


4.1E+15 


l.lE-8 


218.475637 


H2CO 


32,2 — 22,1 


68 


1.8 


8.9 


7.9E+14 


2.2E-9 


226.659543 


CN 


935/1 _ . 13/1 
2 2 2 2 


16 


1.9 


9.6^ 


3.3E+14 


9.0E-10 


226.874764 


CN 


5 5/7/3 1 3 3/5/1 


16 


2.6 


12.3^ 


1.6E+14 


4.6E-10 


658.553275 




6-5 


111 


35 


153 


6.0E+16 


1.6E-7 



^Integrated over channels with emission from 59.4 to 69.4 kms~^. 

^Assuming = 96 K and a nebular mass of 47 M©, equivalent to A^H2 = 3.7 x 10^^ cm~^. 

^Possibly affected by the lack of short baselines and A^x should be considered as a lower 
limit. 
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07'^44"'52f4 52f2 52f0 51 f 8 

RA (J2000) 

Fig. 1. — SMA 654 GHz line- free continuum map {thin lines) toward G240 with the H2 
emission features {gray thick lines; Kumar et al. 2003). The continuum emission is resolved 
into two clumps, with clump A coincident well with a VLA 6 cm point source {cross] Hughes 
& MacLeod 1993). Triangles represent the locations of H2O masers (O'.'l positional error; 
Migenes et al. 1999) while the square indicates an OH maser spot (0^.'4 postional error; 
Caswell 2003). Two pairs of arrows indicate the two outflow axes that have been observed 
in CO (3 — 2) (Hunter 1997): a prominent outflow axis at P. A. = 138° {black arrows) and 
a possible second outflow axis at P. A. = 101° {gray arrows). Contour levels correspond to 
(-2, 2, 4 to 20)x0.13 (la) Jybeam-^ with a beam size of V!5 x O'.'S (P.A. = 6°). 
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Fig. 2. — (a) (6 — 5) velocity-integrated brightness map (contours) from 59.4 to 

69.4 kms"-*^ superposed on the 654 GHz continuum map {gray scale). The C-'^^O emission are 
detected toward both continuum clumps. Symbols follow the use of Fig. [1] with additional 
white crosses indicating the continuum peaks at 654 GHz. Contour levels correspond to 
(-3, -2, 2, 3 to 7)x20 (la) Kkms-^ (6) and (c) C^^O (6 - 5) spectra (histograms) toward 
clump A and B with Gaussian fits {solid curves). 
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Fig. 3. — (a) 223 GHz line- free continuum map. Contour levels correspond to (—4, —2, —2, 
20, 40 to 120) x2 (Icr) mJy beam"^ with a beam size of 4'.'4 x 2'.'4 (P.A. = 1°). White crosses 
indicate the positions of the two 654 GHz continuum peaks. The H2 emission features are 
shown in gray contours (Kumar et al. 2003). (6) through (i) Velocity-integrated, from 59.4 
to 69.4 kms~^, brightness maps of lines in the 223 GHz band (contours) superposed on the 
223 GHz continuum map {gray scale). Contour levels correspond to (—8, —6 to —2, 2, 4 to 
14)x0.8(lcr) Kkms"-^ except those in (e), whose contour levels correspond to (—4, —2, 2, 
6 to 26) X 0.8 Kkms"-^. (j) Position- velocity diagram of a slice through the two continuum 
peaks, where position is the center of clump A. Contour levels correspond to (—4, —2, 2 
to 22)x0.2(la) K. 



